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Abstract 

A four-index tensor is constructed with terms both quadratic in the Riemann 
tensor and Unear in its second derivatives, which has zero divergence for space- 
times with vanishing scalar curvature. This tensor reduces in vacuum to the 
Bel-Robinson tensor. Furthermore, the completely timelike component referred 
to any observer is positive, and zero if and only if the space-time is flat (exclud- 
ing some unphysical space-times). We also show that this tensor is the unique 
that can be constructed with these properties. Such a tensor does not exist for 
general gravitational fields. Finally, we study this tensor in several examples: 
the Friedmann-Lemaitre-Robertson-Walker space-times filled with radiation, 
the plane-fronted gravitational waves, and the Vaidya radiating metric. 

PACS Numbers: 04.20.-q, 04.20.Cv 



I Introduction 

The investigation of conservation laws in general relativity has a long history. From 
its very beginning much of this research was based on pseudo-tensors instead of fully 
covariant methods. The aim of many of these works was to find differential laws which, 
once transformed into integral ones, were interpreted as energy balances in such a way 
that expressions for the energy and momentum densities of the gravitational field could 
be identified. 

*Also at Laboratori de Fi'sica Matematica, Societat Catalana de Fi'sica, lEC, Barcelona. E-mail: 
mangel@fTn .ub.es, cfs@ tpi . uni-j ena.de 
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The covariant approach to this problem is mainly based on the analogy of the 
Bel-Robinson tensor [|| with the energy-momentum tensor of the electromagnetic 
field (there are other approaches based on this analogy, see for instance [^). The 
Bel-Robinson tensor is conserved in vacuum, completely symmetric and traceless. 
Moreover, the completely timelike component referred to any observer (described by 
a timelike unit vector field) is non-negative, and its vanishing implies that the space- 
time is conformally flat (flat in vacuum). This is a desirable positivity property for 
any candidate to gravitational energy density. In spite of these good properties, the 
Bel-Robinson tensor has dimensions of energy density square and this fact makes its 
interpretation somewhat unclear. Nevertheless, it has revealed as a very useful tool 
in many kinds of studies, which has led to some efforts in finding extensions of the 
Bel-Robinson tensor for more general cases than vacuum. Therefore, the question 
arises whether generalizations of the Bel-Robinson tensor exist for space-times not 
necessarily empty. 

The Bel tensor was the first attempt on this problem. It is a tensor whose com- 
pletely timelike component is positive and zero only when the space-time is Minkowski. 
In vacuum, it reduces to the Bel-Robinson tensor, but in the general case it is no longer 
conserved. 

For general space-times, Sachs found a divergence-free tensor that coincides 
with the Bel-Robinson tensor in vacuum. Unfortunately, this tensor does not satisfy 
any positivity property and it is neither completely symmetric nor traceless. 

The systematic treatment of this problem was made by Collinson , who found all 
four-index divergence-free tensors with terms either quadratic in the Riemann tensor 
or linear in its second derivatives. The result found out is that any such tensor can 
be derived from only one tensor, namely T^q'^^ (see (|T7p in the appendix A), whose 
divergence with respect to the first index vanishes. However, there is not enough 
freedom to construct a tensor with its time component positive. 

In this paper we show that, unlike the general case, for space-times with zero 
scalar curvature (i? = 0) it is possible to construct a unique generalization of the 
Bel-Robinson tensor. We begin in Sect. 2 by proving that, when R = 0, there exists 
another conserved tensor which cannot be derived from the Collinson one. In Sect. 3 
we show that, demanding symmetry in the three free indices, there is not any other 
tensor independent from these two. This new tensor allows us to construct (Sect. 4) a 
divergence-free tensor which has the completely timelike component non-negative and 
zero only when the space-time is fiat (excluding some cases that, via Einstein's field 
equations, have an unphysical matter content). This tensor is completely symmetric 
in its three last indices, but it is impossible to get a similar tensor symmetric in all 
their indices. We remark that it is not possible to construct any other tensor with such 
characteristics. 

In order to illustrate this development, we study in Sect. 5 some examples in 
which we can define this tensor, namely: the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) models with a energy-momentum content of (incoherent) radiation {p = f?/3), 
the plane-fronted gravitational waves with parallel rays {pp waves), and the Vaidya 
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radiating space-time. 

Finally, we recall that for purely electromagnetic space-times, and supposing that 
the Einstein field equations hold, Penrose and Rindler also gave a generalization of 
the Bel-Robinson tensor by using spinor methods. This tensor is conserved, completely 
symmetric and traceless in its three last indices. In the appendix B we find its tensorial 
expression and a new (up to our knowledge) positivity property. 

II Deduction of the new conserved tensor for R = 

Otherwise said, throughout this paper we will consider the metric tensor g^is to have 
signature (—,+,+,+). The convention for indices on the Riemann tensor that will be 
used is defined through the Ricci identities: 

(V„V^ - V^V J vx = -R\o.pVa , (1) 

where Va is an arbitrary 1-form. The Ricci tensor and the scalar curvature are defined 
as usual: Rai3 = R"aap ^'^d R = R". We also recall the Riemann symmetries and the 
first and second Bianchi identities: 

RapXfi = R[al3][\fjt\ = RxfiaP i 
R[al3X]fi = , 

"^[uRaflXfi = . (2) 

The procedure we are going to use to find the conserved tensor starts from the 
expression for the divergence of the Bel tensor |[l[]: 

where the Bel tensor T""!^^^ and J°^-^ are defined as follows: 

rj^c^pxn ^ 1 [r^p^'^rI^^ p^ + *R*''P^''*R*l^p P^ + *R"P^''*Rf^p P^ + R^'^P^" R^f^p , 
JA/./3 ^ y>^Rf^P - V'^i?^^ = V^i?^^^'" , (4) 

being "*" the usual dual operator acting over any pair of antisymmetric indices: 

*Ral3Xp, = \'<]al3upR 'a/x 5 R*apXp. = ^VXpapRa/S ^ ! 
*R*al3Xf, = iVa/SySVXfiapR"'^''^ , 

and rjaisxp is the canonical volume 4-form. Our purpose now is to work out the right 
hand side (r.h.s.) of equation @ in order to convert it into a global divergence. To 
that end, we will repeatedly integrate by parts and make use of equations (^0). 
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To begin with, the last term in equation can be easily transformed into a 
divergence, by means of the Ricci and Bianchi identities (|l],|D: 

-]^g^^RP^^^r^f = -^Vv^VpJ""^ . (5) 

Next we expand the leading terms of the r.h.s. of by using the definition of J"^^ 
and integrating by parts: 

R'^p ^^V^/?"" - V„ [r^'pR^^ + R^PV^R^p ^" + [A ^ — ^ /i] . (6) 

The first term in the previous expression can be rewritten by means of the Ricci 
identities (|l|) as follows: 



^ r , 



V" [r^'^\R'"^) + R'^^V'Ri^^] + R^V^R'"^ . (7) 



Thus, we have converted the r.h.s. of equation (|^) into a divergence plus the following 
terms: 

= i?"^ (V^R^^ - V^R^) + R^'^V^R^ + R^P (VpRf^^ - V^Rp) + [A < — > fi] . (8) 

Now, taking into account the contracted Bianchi identities {V p,R^^ = ^V^R), these 
terms can be transformed into the following expression: 

-V^ {r^^R""^) + (R^^R"^) + (r^^R"^) - 2R^'' (V^R^ + V^R^) + 

[r^^R^" + R^^'R^P + 2i?^"i?^^) - i (^R^^V^R + R^^V^R + 2R^''V^r) . (9) 

The last three terms of this expression vanish when R is constant, so we are finally 
left with —2R^" {v^R^ + V'^-R^^ . Nevertheless, notice that our final purpose is to 

find a conserved tensor j"""^-^^ whose completely timelike component referred to an 
observer u, T"°'^^'^UaUf3UxUp_, is positive, which means that we are only interested in 
the symmetric part. Therefore, without lost of generality, we can symmetrize the 
whole expression and, as a consequence, the remaining terms transform themselves 
into a divergence: 

_4^'^(/3v^i?^) = -2V(^ [r^R^^'') . (10) 

So, we have finally achieved a conserved tensor if the scalar curvature vanishes (in 
fact, if it is constant). Collecting all the previous terms ([5|-|T0|) we get the final result: 
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where we have defined 



- 2g'^^^V^V^R^^'' - V „V R"^^ g^'''^ + V R""^^ g^^'^ . (11) 

It is a matter of checking that this tensor cannot be obtained from Colhnson tensor 
TiQ^^ when it is restricted to the case i? = 0. Now, we are left with the question of its 
uniqueness. 

Ill Uniqueness 

In this section we will prove that, in the case we are concerned with (i? = 0), it does not 
exist any other conserved tensor, symmetric in its three last indices, independent from 
rpna/sxtj. ^j^g Collinson tensor, T^^^^. The reasoning is the following. Suppose you 
are given a tensor y'""^-^^ which is conserved when R = 0. The divergence computed 
in the general case will be a combination of the following type (taking into account 
symmetries and unit dimensions): 

V«T'"°^^^ = aR^f^^V^^R + bRV^f^R^^^ + cR'^^'^V^Rg^"^ + dRg^f^^V^^R + 
+eV^^V^V^^R + fVV^'^Rg^'''^ + hV'V^^V^Rg^^^ + iV^^V'"' V^i?^^^) , 

a, b, c, d, e, /, h and i being constants. This can be immediately cast in the following 
form: 

V^T'""^^^ = {a- b)R^''^W^R + V„r"^^^ , 

where r"^'^^ stands for: 

^a/3AM = bRg''^'^ R^i"^ + c (^i?/?"(^c/^^) - i^"('3c/^'')i?2^ + c/i^(^^c/'^)"/22 ^ 

+e^"(/3yAy^.)^ + f g'^^'V aV^^ Rg^^^ + hg'^'V^^V ^Rg^'''^ + ig'^^^V^^'^V aRg^^'^ , 

so clearly it is a tensor that vanishes when R does. That is, if a tensor of the kind 
we are considering is divergence-free when i? = 0, in the general case its divergence 
should be a multiple of R^^^V^^'R plus the divergence of a tensor of the type t""^^^. If 
we had two such tensors, a suitable combination of them removing the term R^^^V^'^R 
would give a conserved tensor for the general case and, therefore, due to the Collinson 
result 0, it could be constructed from Tiq'^^. Given that the tensor r"^^'^ vanishes 
when R = 0, the three tensors would not be independent in that case. 

On the other hand, this reasoning shows that, from the very beginning, we were 
able to know that in the R = case at most one more conserved tensor could exist 
apart from Collinson's one, as finally has been the case. 
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IV Positivity 



As it has been pointed out above, in the general case all the conserved tensors can 
be constructed from Tiq^^. This construction is based in two procedures. First, it is 
clear that if we perform any permutation on the three last indices we will still have a 
conserved tensor. Second, by taking the two traces T^i^'p, ^icT^p and multiplying them 
by g^^ we obtain new conserved tensors. Actually, there is not any other two-index 
divergence-free tensor independent from them. These two tensors can be taken to be 
(as usually obtained by hamiltonian differentiation): 

tf = 2V''V^R-2g'"^Vf,V^R+^g"'^R'^ -2RR''^ , 

tf = 2V^V"i?"^ - V,V"i?°^ - 2i?°"i?^ + ^R^pR'^Pg'"^ - ^^"^V^V^i? . 

This is all the freedom we have in the general case. 

In the R = case, there exists only one two-index conserved tensor, which is that 
obtained from tg'^. On the other hand, as we consider tensors which are symmetric 
in the three last indices, from the CoUinson tensor we will only have one independent 
tensor (apart from the two-index tensor), namely Tiq^^^'' or, equivalently, the Sachs 
tensor T'°'^^^ which is symmetric in its three last indices. 

Therefore, to construct the super-energy tensor in the R = case, we are led with 
three tensors: i) the Sachs tensor T'"^'^'^ (restricted to the R = case): 

j.,af3Xi, ^ ya(/3AM) ^ 2i?°"i?^^^^^) - - -R^pR'^fg^^^g^^'^ + 

3 2 

+ 2V^^ V^/?^^" - - V^^V'"! i?^'') - - V" V^^ R^^^ + 
3 3 

+ -V,V"i?(^^^^)" - 2^"(^V^V^i?''^" - V^V'^R^^'^g^"^ , 
3 



J the two-index tensor t"^: 



= 2V^V"i?"^ - V^VR"'^ - 2R'"'R1 + ]^R„pR''Pg'''^ , 

and in) the tensor T"""^^^ previously found in (|TT]). 

Now, we have to combine these three tensors in such a way that any observer 
measures a positive quantity. Moreover, we would like that this completely timelike 
component vanishes if and only if the space-time is flat. 

First of all, we have to take into account that terms made of derivatives of the 
Ricci tensor do not have a definite sign, so it would be necessary to eliminate their 
contributions. This aim can only be achieved by means of the following combination: 



~ 9 V J 2 
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which more exphcitly reads: 

- 2^''(^V^V^i?^)'" + 2V^V'"i?(^^^^)" - 2V^V''i?°(^^^^) . (12) 

Since we have aheady exhausted all the freedom, we finally examine the completely 
timelike component referred to any timelike unit vector u: 

A{u) = A^p^^u''u^u\^ = T^py^pu''u^u\^' + ]^R,pR''P + 

+ 2 (Ro^pu'^u")' + 3 (R'^'R^) u^Ufi - iVaS/'^R'^) u^up . (13) 

To check the positivity of A {u) it is convenient to write out the last term of the previous 
expression in the following form: 

V,V"i?"^ = C^fR"" + 2i?""i?^ - ^g'^^R^pR'"' . 

At this point, we introduce four spatial tensors, namely Eapiu), Hapiu), Ma/siu) and 
Na/siu), that (together with R) wholly characterize the Riemann tensor |^. Their 
definitions, properties and some useful formulae are given in the appendix A. 
Introducing the previous definitions in (13) we obtain, after some calculations: 

A (u) = {E^p - M^p) {E"" - M'^P) + H^pH^" + 3iV^,Ar'^'' + , (14) 

which is a sum of square terms (all the tensors appearing here are spatial). Therefore 
it is manifestly positive and its vanishing implies: 

HafB = , Nap = , 

The previous expressions lead to Rapu"' = (see appendix A). This condition, when 
considering the Einstein field equations, immediately drives to an unphysical energy- 
momentum tensor. Hence, if we eliminate these unphysical space-times (for instance, 
adding any energy condition), the vanishing of A{u) finally implies the Minkowski 
space-time. 



V Some examples 

In this section we are going to study the tensor A'^'^^^' (|g) in some space-times with 
vanishing scalar curvature. In particular, we are going to consider the following ex- 
amples: (i) the radiation FLRW cosmological models, (ii) the pp waves, and (iii) the 
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Vaidya radiating metric. In examples (i) and (ii) we will give the expression for A"'^^^ 
and its completely timelike component A{u) ([l^) for an arbitrary observer u. In the 
example (iii), for the sake of brevity we will give only the expression of A{u), also for 
an arbitrary observer. 

In the first example we study the case of the FLRW models (see for instance pf) 
with vanishing scalar curvature, the radiation models, whose energy-momentum tensor 
(of perfect-fluid type) is given by the following expression (throughout this section we 
will use units in which SttG = c = 1) 

where U is the fluid velocity {U"Ua = —1), Q the energy density, p the pressure, and 
hap = Qap + UJJp the orthogonal projector to the fluid velocity. The line element of 
these conformally-flat models can be written as 

ds^ = -dt^ + a\t) {dx^ + S2(e, x) {dO'^ + sin^ ^V)} , 

where S(e, x) is given by 

sin X if e = 1 , 
E(e,x) = { X if e = , 
sinh X if e = — 1 . 

The fluid velocity U, the scale factor a(t) and the energy density g(t) are 

d SA^ 
U = —, a\t) = {t- to) [2A - e{t - Q] , g{t) 



respectively, and A and to are arbitrary constants. 

After some straightforward calculations, and using the special properties of these 
space-times, we arrive at the following expression for A""^^^ 

3 L 3 3 J 

As we can see, it is proportional to the energy density squared. We can also check that 
it is indeed divergence-free. Now, let us compute the completely timelike component 
(0) of this tensor with respect to an arbitrary observer u. To that end, we decompose 
u in the next way 

u = -f{U + v), w"f/„ = 0, v''va = v^>0, 7 = (1 - i;2)-i/2, 

where the case v = corresponds to an observer comoving with the fluid {u = U). 
Then, from (|1^,|15|) we find that A{u) is given by 

A{u) = ^g'j^{l + lv' + lv^ 



3" ' I 3 3 
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That is, it is function of q and v only. Moreover, it increases monotonically as v 
increases and its minimum corresponds to the case w = 0, in which the observer is 
comoving with the fluid. 

Now, we are going to study the tensor A""^^^ in the case of the pp waves space- 
times. The corresponding line element can be written in null coordinates {m, f , C, C} 
follows (see |^ for details) 

ds^ = -2dudv + 2dCdC, - 2Hdu^ , 

where H is an arbitrary function which does not depend on f [if = H{u, Using 
the following Newman-Penrose basis {£, fc, m, m} 

-a - (9 (9 d 

I = — , k = — -n t;— , = TT 1 

ov ou ov oC, 



the Ricci and self-dual Weyl tensors are 

Rap = 2$£q,£/3 , Caf3Xn = CafiX^i + i Caf3Xn= 2\I/4V^/3Va^ , (16) 

respectively. Where the quantities \I^4 and Va/s are given by 



From (|Ig) we can see that the energy-momentum content can correspond with vacuum, 
Einstein-Maxwell or pure radiation fields. Moreover, the Petrov type is N, being i the 
repeated principal direction of the Weyl tensor, which in fact is a constant vector field 
(V«£^ = 0). 

After some calculations, we have found that the Bel tensor T""^^^ and our tensor 

From this expression, the completely timelike component is given by 

A{u) = 4(3$2 + ^4^'4)(rM„)^ . 
Then, the vanishing of A{u) implies the Minkowski space-time. 

Finally, we are going to consider the Vaidya radiating space-time (see for instance 
0). For the sake of brevity we only give here the completely timelike component ([T^) 
of the tensor A°'^^^^. The line element of this spherically symmetric metric can be 
written as follows 

ds"^ = -2F'^{u, v)dudv + r^{u,v){de'^ + sin^ dd^"^) , 
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where 



dr 



dr 



F\U,V) = fiu)—, T^ = 7T/(n 



'2m{u) 



dv 



du 



- 1 



Here, m{u) is the invariantly defined mass function. As is well-known, the Petrov type 
of this metric is D. Then, taking the following Newman-Penrose adapted basis 



-ld_ 
'Ydv 



k 



-ld_ 
'Ydu 



m 



d id 
+ 



where ^ and k are aligned with the principal directions of the Weyl tensor, the Ricci 
tensor is given by 



R. 



a/3 



m . 



And the only non-zero component of the Weyl tensor in this basis (see |^) is 

m{u) 

$2 : 



r2('U, v) 



which in this case is real. 

In terms of these quantities, we have found the following expression for A{u) 



A{u) 



^2 + 



+ 



36(rM„)2(A;%)2 - 18(rM„)(A;%) + ij + im^{tuc 
and again, taking into account that 

2(rM„)(fc%) > 1 =^ 36(ru„)2(fc%)2 - 18(ru„)(A;%) + 1 > 1 , 

A{u) vanishes if and only if the space-time is the Minkowski space-time. When we re- 
strict ourselves to observers lying on the 2-planes generated by the principal directions 
[2(£"Uq,)(/c'^m^) = 1], which are precisely the observers that minimize A[u), the result 
is 



A{u) 



^2 + $(rM«) + 4^2 + lo^irw 



Acknowledgments 

M.A.G.B. gratefully acknowledges the Comissionat per a Universitats i Recerca de la 
Generalitat de Catalunya for financial support. C.F.S. gratefully acknowledges financial 
support in the form of a fellowship from the Alexander von Humboldt Foundation. 



10 



A Useful definitions 

We first write down the Collinson tensor |^: 

Tfcf = 6Qf + gf ^^'^^ + gf ^^^^ , (17) 

where 

Qf^" = -Sg^'^RJ^RP^P'' + 8R^^^pR'^''"P + SR^R^''^" + mf^pR^^'^P + 
+2i?^i?^°"^ + 2R/^pR^''''P + 3^°^i?^i?^" - g''^R^p\RP^^'' , 

which is divergence-free in the index a and whose only symmetry on the indices /5, A 
and /X is Tfo^l'^l - TiXt^'^ = 0. 

Next, in order to introduce other useful definitions, recall the well-known decom- 
position of the Riemann tensor into its irreducible parts under the full Lorentz group: 

where Ca/3Xfi is the Weyl tensor and 

Ea/BXp. = - {RaX9f3fi ~ R-a^iQpX + RpfiQaX — R^XQai^j i R-a/B = Ra/B ~ -^RQoP i 
GafBXn = — {gaXgpiJL " dafig^x) , 

being R = R^ the scalar curvature. 

The electric and magnetic parts of the Weyl tensor associated to a timelike vector 
field u are: 

Eaxiu) = CapXtiU'^Up , Ho,x{u) = - C afiXi, U^U^ . 

These tensors are spatial (orthogonal to u), symmetric and traceless, and they fully 
determine the Weyl tensor. 

We can proceed analogously with the tensor E^ibx^i and define (see [0 for a more 
detailed study of these matters): 

We give here some of their properties, since they are less known than the electric and 
magnetic parts of the Weyl tensor: 

M,a = Ma„, M,am^ = 0, Mli = R^,uV, 
Nax = -Nxa, N^xu^ = 0, iV/^ = 0. 
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The tensor Max has 6 independent components, while Nax has only 3. Actually, they 
completely characterize the traceless Ricci tensor: 

Rap = -2Map - 4 * Ar,(„M"M/3) + M; {gap + lUaUp) . 

From the previous definitions, it is clear that 

NaP = R<.pU''u^ = R = ^ =^ Rap = -2Map, 

and this implies, in particular, that Rapu^ = 0. 

Let us finally give some formulae which is useful for the derivations of some expres- 
sions in Sect. 4: 

T^^^^'UaUpUxU^ = EapE"^ + H^^Hap + M,^M°^ + NapN'^^ + ^ , 

(i?„X)(^X) = 2iv.,iv-'' - (m;)% 



B Electromagnetic case 

For space-times with an electromagnetic energy-momentum content, Penrose and 
Rindler gave the following modification of the Bel-Robinson tensor (see for 
the spinor notations and conventions): 

tapXfj. = '^ABCD'^A'B'C'D' — 

-2^^/ CD'^AB^ C>D^A>B> + D{A'^{AB^ C)\D'^ B'C) ' (18) 

where abcd and ^ab are the Weyl and the electromagnetic spinors respectively, and 
7 is the gravitational constant. This tensor is symmetric and traceless in the first three 
indices and has zero covariant derivative with respect to the last one, provided that 
the Einstein field equations hold: 

^^tapxy. = , 
tapXfi = t(apX)fj. , taXfj. = • 

It is important to notice that the second part of this tensor is formed with (VF)^ terms, 
so it cannot be expressed by means of the Ricci tensor and its derivatives. Therefore 
it is independent from the tensors considered above. 



The tensorial expression for (18) is rather involved but, in this case, it is useful in 
order to proof a positivity property. It should be noted here that in the final result 
we have returned to the initial signature metric (—,+,+,+), and units such that 
8717 = c = 1. With these conventions: 

tapXfi = -^apXti + 2V(Q,F|^|_aVA)-F'^ + 2V(aF|^|^V'^Fx)/, - 

-\^[aF\„p\^pF''pgx), + V.F,(„V''F^^^7a)m - Qi^ap^'' F'^^V pF,^ , (19) 
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where T^px^ is the Bel-Robinson tensor and Fq,^ is the electromagnetic tensor. From 
this expression it is easily seen that ([T9| ) satisfies the following positivity property: 



t^px,u uuu - U ^ I ^ u'^W^Fx, = 0. 

To prove this, let us introduce the orthogonal projector to u, h^p = Qap + Then, 
we compute ta/3Xfj.u°'u'^u'^u'^: 



= \%px,u'^u^u\'' + {h^'^F^pu") [hx^F^^u^) + 
+ 1 [h^-h^Pp^p] {hx^KpP^^^) > , 



2 

and the equality holds only when Capxfi = and 

h'^^h^^Px, = , h'^'^P^pu^ = 0^p^p 
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